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Abstract. When Einstein's equations for an asymptotically flat, vacuum spacetime are 
reexpressed in terms of an appropriate conformal metric that is regular at (future) null 
infinity, they develop apparently singular terms in the associated conformal factor and 
thus appear to be ill-behaved at this (exterior) boundary. In this article however we show, 
through an enforcement of the Hamiltonian and momentum constraints to the needed 
order in a Taylor expansion, that these apparently singular terms are not only regular 
at the boundary but can in fact be explicitly evaluated there in terms of conformally 
regular geometric data. Though we employ a rather rigidly constrained and gauge fixed 
formulation of the field equations, we discuss the extent to which we expect our results 
to have a more 'universal' significance and, in particular, to be applicable, after minor 
modifications, to alternative formulations. 



1 



I. Introduction 



The natural outer boundary for gravitational radiation problems in asymptotically 
flat spacetimes is future null inflnity (also referred to as X+ or 'Scri') [1,2,3,4]. The 
feasibility of using hyperboloidal slicings of spacetime and putting the outer boundary of 
numerical calculations at Scri has long been recognized, thanks to the elegant work of H. 
Friedrich and his associates who have numerically implemented Friedrich's conformally 
regular field equations for this purpose [5,6]. In spite of this success, most numerical 
relativists attempting to solve the black hole collisions problem prefer to use some direct 
formulation of the Einstein equations themselves (which, though implied by Friedrich's 
conformally regular system, are not explicitly included within it). Since, however, 
Einstein's equations are not conformally regular in any obvious sense, it has heretofore 
not seemed feasible to put the outer boundary of such, more conventional numerical 
calculations at Scri. 

In this paper however, we propose a specific (constrained and gauge fixed) formu- 
lation of the Einstein equations for which it does seem feasible to put the exterior 
boundary at future null infinity. In particular, we show that when the usual Hamilton- 
ian and momentum constraints are enforced (at least to a suitable order in a Taylor 
expansion about Scri) and when the well known geometrical condition that Scri should 
be shear-free is taken into account, then all of the apparently singular terms in the Ein- 
stein evolution equations (that appear when the latter are reexpressed in terms of the 
conformal metric) are not only actually regular at this boundary, but also can in fact 
be explicitly evaluated. We further show that the conditions at Scri that were needed 
to establish this regularity (i.e., the vanishing of the shear and satisfaction of the con- 
straints to the requisite order) are preserved by the boundary forms of the evolution 
equations so derived. 

Our particular formulation of the field equations was motivated, in part, by the 
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proof, given in Ref. [7], that the usual vacuum Einstein equations expressed in CM- 
CSH (constant-mean-curvature-spatial-harmonic) gauge satisfy a well-posed Cauchy 
evolution theorem in the 'cosmological' setting of spatially compact spacetimes. The 
constant mean curvature gauge condition serves to partially decouple the constraints 
and to define a convenient slicing for spacetime whereas the spatial harmonic condition 
effectively reduces the evolution equations for the spatial metric to hyperbolic form. 
The remaining spacetime metric components (i.e., the lapse function and shift vector 
field) are determined by elliptic equations that enforce the continuation of the gauge 
conditions. This theorem, on the other hand, does not treat the decomposition of the 
spacetime metric into the product of a conformal metric and a conformal factor nor 
does it address constructs appropriate to asymptotically flat spacetimes such as the 
attachment of a conformal boundary. Thus it does not literally apply to the formu- 
lation proposed here. Nevertheless its existence suggested the feasibility of adopting 
analogous gauge conditions (wherein the physical mean curvature is held constant but 
the conformal metric is instead subjected to the spatial harmonic condition) for the 
present problem. 

In spite of the built-in rigidity of our particular setup, we are confident that the 
central results of the present paper have a more 'universal' significance and do not 
require, for their application, a strict adherence to our specific formulation of the field 
equations. Rather we believe that our main results can be applied to a wide variety 
of alternative formulations that employ, for example, different gauge conditions or 
allow free (i.e., unconstrained) evolutions provided suitable attention is paid to the 
necessary regularity conditions at Scri. We shall return to and clarify these remarks 
about 'universality' later in the presentation. 

Though a full numerical implementation of our ideas has not yet been carried out, 

one of us (O.R.) is planning to modify his existing vacuum, axisymmetric evolution 
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code [8] so as to fix the outer boundary at Scri using a variant of the ideas described 
herein. Furthermore, L. Buchman and H. Pfeiffer have estabhshed the feasibihty of 
numericaUy solving the constraint equations out to Scri by exphcitly computing initial 
data sets for multiple, boosted spinning black holes [9] . In addition the linearization of 
our proposed formulation, about a conformally compactified Schwarzschild background, 
has been carried out and partially analyzed with a view towards developing code tests 
for the nonlinear calculations [10]. 

II. Constraint and Evolution Equations 

Our constrained evolution scheme involves solving elliptic equations for the confor- 
mal factor, the conformal lapse and physical shift and for the unphysical (i.e., confor- 
mal) mean curvature. The conformal spatial metric and a momentum tensor density 
(the trace-free part of the physical ADM momentum) are evolved. The constant time 
slices are CMC slices which extend outward to future null infinity which is fixed to 
coincide with a coordinate cylinder in the conformal spacetime metric through the 
imposition of suitable boundary conditions upon the conformal factor, the conformal 
lapse and the physical shift. 

The physical ADM variables are {^y,7r*-'} and {N,X^}, where gij is the spatial 
metric, tt*-' is a momentum tensor density (i.e., tv^^ / Hg, with Hg — ^/det~g^, is a 
tensor) defined in terms of the extrinsic curvature Kmn via 

= i^gig^'g'^'' - g'^gnKmm (2.1) 

N is the lapse and X* is the shift. In our conventions JCrnn is defined via gmn,t = 
—2NKmn + X'^gmn,i + ^^mgin + ^%gmi- We usc Latin (Greek) indices to denote 
spatial (spacetime) indices. The spacetime metric has the line element 

= -N^dt^ + gij{dx' + X'dt){dx^ + X^dt) = ^^^Qf^^dx^dx" (2.2) 
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and we write this alternatively as 

2 _ {-N'^dt^ + -fij{dx' + X'dt){dx^ + XHt)} _ (^^T^^.dx^da;'^ 
ds - - , (2.3) 

where il is the conformal factor, 7^-, the conformal spatial metric and N the conformal 
lapse. Notice that X* serves as both physical and conformal shift but that gij , , A'" 
and N are related by 

1 N 
!»J = N = -. (2.4) 

Q is greater than zero in the interior but approaches zero at X+. 
CMC slicing is defined by 

= -2K = constant (2.5) 

where tr^Tr = gijir^^ and K = —g^jK"^^ = —(mean curvature). The slightly odd sign 
convention chosen for the definition of K has been made so that this constant will be 
positive. The traceless part, n^^'^^ , of the ADM momentum is defined by 

„tri7 . „i7 ^ „ii „ ^rnn /n «\ 

TT J -—T^ J - -g Jg^^-K . (2.6) 

When K is Q. constant the momentum constraint can be written 

tr ij ■ 



(2.7) 



where is the covariant derivative with respect to 'jij. Defining, as usual, the mixed 



physical components of tt*'' via irf^-' = gim'^^^^'' the momentum constraint can also be 



tr j 

jiiipwiiciito KjL n via /'' 

written as 



V,-7rf^' = 0; nr' = 0. (2.8) 

Standard methods are available for solving this equation through the use of orthogonal 
decompositions. 



The Hamiltonian constraint (also known as the Lichnerowitz equation) has the form 

- 4i](v^ VjQ) + GY^n^i^j - R{-f)n'^ (2.9) 

where V-' Vj = 7*-' VjVj is the Laplacian with respect to 7^- , -^(7) is the scalar curvature 
of 7jj and //^ = i/det 7mn- This elliptic equation for Q degenerates at X+ where the 
conformal factor vanishes. 

From the ADM evolution equations, one can easily derive (by taking the trace of 
the Qij^t equation) the following evolution equation for fi, 

6[n,t-x'n,i] -nv = -2KN, (2.10) 

where F is defined by 

r:=7'""7mn,t-2V/X'. (2.11) 

To fix the decomposition of gij into a conformal metric and a conformal factor, we need 
to impose a normalization condition upon 7^-. A mathematically appealing choice is 
to exploit Yamabe's theorem [11] and demand that -R(7) be a (spacetime) constant. 
The corresponding requirement that dtR{'y) = then leads to an elliptic equation for 
F given by 

^ Y^ (v. v,r) + i ^(7) F = V, V, 

where Rij (7) is the Ricci tensor of 7^- and, as defined above, N is the conformal lapse 
function {N := NQ). Note also, from Eq. (2.11), that the quantity -F/(2iV) is the 
unphysical mean curvature. 

N is determined by solving the elliptic equation 

Q TV T 

^ ^ 2b ^2.13) 

^-^W) + I ^ lilljmT^ TT 
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trij 



-Rijii) 



2N 



IT 



trij 



(2.12) 



which enforces the condition dtK = (assuming that K is spatially constant). By com- 
bining Eqs. (2.9) and (2.13) in a straightforward way, one can derive a less degenerate 
form for the conformal lapse equation which has only a single power of VL multiplying 
the Laplacian of N. 

For the physical shift vector X*, we propose to determine it so as to preserve the 
spatial harmonic gauge condition defined by 

^':=7^^(f^,-(7)-f^,(7)) = (2.14) 

where 7 is a fixed (i.e., time-independent) reference metric (for example 7ij = 7^ |t=o is 
a possible choice). The Christofi^el symbols of 7 and 7 are denoted by f ^^(7) and f ^^ (7), 
respectively. Equation (2.14) corresponds to the demand that the identity map from 
(M, 7ij) to (M, 7y) be harmonic. The shift equation results from requiring dtV^ — 
and is given by 

^ = -\rv^-\i^'r,i - ^TT*"^- {t%{^)-v%{^)) (2.15) 



^7 



- V* 



(V^X^- + V^X^) (f^,-(7)-f^,(7)). 



An extensive mathematical study of the use of constant-mean-curvature-spatial- 
harmonic (or CMCSH for brevity) gauge conditions in a 'cosmological' (i.e. spatially 
compact) setting was made by the authors of Ref. [7] who proved a well-posedness 
theorem for the vacuum Einstein equations (in arbitrary spacetime dimension) in this 
gauge. Their theorem did not deal with the conformal decomposition of the spacetime 
metric or with the presence of future null infinity in the conformal geometry and so, 
strictly speaking, is not applicable to the problem dealt with here. Nevertheless, their 
formulation provided some of the motivation for our setup and might conceivably pro- 
vide the model for the well-posedness theorem that one would eventually like to prove 
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for the equations we are studying. It is worth mentioning here that the main reason for 
the choice of the spatial harmonic gauge condition made in Ref. [7] was the fact that it 
nullifies terms in the Ricci tensor of gij that, if present, would disturb the hyperbolic 
character of the equations of motion for this metric. The strategy is rather similar to 
that motivating the use of spacetime harmonic coordinates in other formulations (to 
nullify corresponding terms in the spacetime Ricci tensor) but leaves open the possibil- 
ity for determining the lapse and shift through the solution of elliptic equations instead 
of hyperbolic ones. Of course the resulting evolution system is now hyperbolic/elliptic 
rather than purely hyperbolic. In our formulation the conformal metric 7^- (normalized 
by the condition R{'j) = constant and gauge fixed by the spatial harmonic conditions 
V'' = 0) represents the two propagating degrees of freedom of the gravitational field. 
At X+, where Q vanishes, Eq. (2.9) forces the gradient of Q to satisfy 

|x+= ) (2.16) 

whereas Eq. (2.10) yields 

X'^i |x+= ^ KN \j+ . (2.17) 
The shift field at must therefore take the form 

X' |x+= ^ NY'nj I1+ \x+ (2.18) 

where the vector field Z* \x+ is purely tangential to (i.e., satisfies 
Z* \x+ |x+= 0). Note that the above then implies 

^ijX'X^ (iV^ + ^ijZ'Z^) |x+ . (2.19) 

If, however, we wish to have the time coordinate vector field ^ tangent to the null 
generators of (i.e., to be such that ^ is null at X+ in the conformal geometry) then 



we would need to have 



^ijX'X^ U+=iV2 \x+ (2.20) 
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satisfied there. In other words, the boundary condition |x+= corresponds to 
choosing ^ \j+ to be null, rather than spacelike, at the boundary. For the 'normal' 
component of X* at X"*" however, one has no flexibility since the above choice is forced 
by the requirement that X'^ coincide with a fixed cylinder in the conformal geometry. 
Finally, the evolution equations for {'^ij, tt*''*-'} are given by 

dtlij = ^ l^l Ijm Tt'"^'"^ + \ 7^i T + Hi Ijm (v'X^ + V^X') , (2.21) 

and 

IN 

d—trii (v'm-.tvii\ -yi ^trmj vj -.trim _trim_tr97., 

-^§7r*"^i^ (2.22) 
+ //^ (v'V^N - i 7"^"V^Vniv) -fiyN (^^^^ (7) - i ^(7)^ 

-^^-""[-^ — s^'^-^r 

The main content of this paper involves analyzing the apparently singular, 0-dependent 

terms in Eq. (2.22). We shall, in fact, derive explicitly regular forms for the evolution 

equations at X+ and show how these can be used to propagate geometric data along 

this conformal boundary. 

In effect, our formulation is maximally elliptic in that we envision solving elliptic 

equations for the Hamiltonian and momentum constraints, Eqs. (2.9) and (2.7), the 

conformal lapse and physical shift, Eqs. (2.13) and (2.15) as well as Eq. (2.12) for the 

unphysical mean curvature. This could only be numerically practical through the use 

of elliptic solvers that operate at the level of 'linear complexity' but fortunately several 

such systems are currently available. However, as we shall emphasize throughout this 

paper, our main conclusions should hold without the need for a strict adherence to the 

constrained evolution program outlined above but should instead be applicable, with 

minor modifications, to a wide variety of alternative schemes. 
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III. Regularity at Future Null Infinity 



The constraint and evolution equations presented in section II, as well as the elliptic 
equation for the conformal lapse function N, are formally singular at Scri where the 
conformal factor vanishes. In this section, we analyze the behavior of the corresponding 
fields in a neighborhood of this boundary and derive a set of regularity conditions that 
sufficiently differentiable solutions must satisfy at Scri. Later we shall establish the 
consistency of the regularity conditions by showing that they are preserved under time 
evolution. This last step will necessitate explicit evaluation of the evolution equations, 
including their apparently singular terms, at the conformal boundary. 

Let M be a three-dimensional spacelike slice with conformal boundary dM S'^ on 
X'^. We choose coordinates {x*} = {{x^,x°') | a = 2, 3} = {(r, 9, (p)} for the Riemannian 
manifold (with boundary) (M, jij) such that = r is a 'radial' coordinate satisfying 
r < r+ = constant on M with r = r+ corresponding to the boundary 2-sphere dM at 
Scri. The indices a,h,. . . are restricted to range over 2 and 3 so that {a;"} = {{6, ip)} 
are 'angular' coordinates for the r = constant surfaces. Aside from requiring that dM 
coincide with the coordinate sphere defined by r = r+, we leave, for the moment, the 
choice of coordinates on M arbitrary. 

We suppose that each of the relevant fields can be expanded in a finite Taylor series 
(with remainder) about the boundary at r = r+. Thus for each u e {fl, 7^, tt*''*-', TV} 
there is an integer I > (depending on the choice of u) such that u can be expressed 
as 

u{x') ^ uoix") + ui{x''){r - r+) (3.1) 
+ ^«2(a;")(r - r+f + ... + ^ui{x''){r - r+Y 
+ rm(a;''), 

with Uk := lim (d^u) and remainder rm(a;*) = o((r — r+Y), on some interior neigh- 

r/'r+ 
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borhood of dM. In the following, it will be convenient to utilize the symbol = to 
denote equality at dM. By definition the conformal factor thus satisfies Q = so that, 
in particular, ilo = in the corresponding Taylor expansion. A remarkable feature 
of the degenerating elliptic equations that we have to deal with is that they permit 
one to explicitly compute more detailed asymptotic information (in the form of Taylor 
expansions) about the corresponding solutions than would be possible in the case of 
non-degenerate equations. In particular, we shall be able to evaluate the first three 
radial derivatives of Q, the first two such derivatives of and the first radial deriva- 
tive of the TT^''^* components of tt*^*-' 'universally' at Scri (i.e., expressible in terms of 
data there without reference to the actual global solutions). Remarkably these partic- 
ular derivatives are precisely what is needed to then evaluate the evolution equations 
(including their apparently singular terms) at Scri and to verify that they imply the 
preservation of the associated regularity conditions. 

A rigorous treatment of the constraint equations on 'hyperboloidal' initial data sur- 
faces intersecting X+ has already been given by Andersson, Chrusciel and Friedrich in 
an important series of papers from the early 90's [12, 13]. In particular they derived 
the regularity conditions needed for a differentiable Scri in terms of ADM Cauchy data. 
Our main contribution here is to carry this analysis a step further and show how one 
can evaluate the (apparently singular) evolution equations at Scri and use them to 
verify preservation of the regularity conditions within the framework of our particu- 
lar gauge fixed evolutionary formulation. Though we use this framework in order to 
have a complete, coherent system for calculations, we do not believe, as emphasized 
in the introduction, that our main conclusions hinge crucially upon its specific form 
but rather that they should apply equally well to a variety of other formulations of the 
field equations. Though the arguments sketched above yield expressions for the first 

radial derivative of only the ri components of n^^^^ at Scri we shall present, in section V 
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below, an alternative method for evaluating the apparently singular terms in the tt*^*'' 
evolution equations that will finally allow us to compute the first radial derivatives of 
the angular components, tt*^"'', at Scri as well. 

As in section II, we assume that the Riemannian manifold (M, 7^) corresponds to 
a CMC slice in the physical spacetime and take the mean curvature of the latter to be 
a negative constant (written as before as —K,K > 0, constant). It is convenient to 
reexpress the conformal metric relative to the chosen coordinates {x^} = {(r, 6, (p)} 
in 2 + 1 dimensional (Riemannian) ADM form by setting 

dl^ = -fijdx'dx^ (3.2) 
= n^dr^ + habidx" + Y''dr){dx^ + Y^dr). 

The induced metric on an r = constant surface thus has the line element 

da^ = habdx^dx^ I , , (3.3) 

lr=constant ^ 

and we let the symbols | o or ' Va(/i) signify covariant differentiation with respect to 
this metric. 

The unit outward pointing normal field to an r = constant surface in {M,jij) is 
given in coordinates by 

1 -¥"- 

= ) (3.4) 

n n 

or, equivalently, in covariant form, by 

(i.O = (n,0,0). (3.5) 

The second fundamental form Xab, induced by 7^ on such r = constant surfaces is 
defined via 

Aab = -VfeJ^a = -l^a;b (3.6) 
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(where, as before, Vj or ; z signifies covariant differentiation witli respect to 7^ ). Writ- 
ten out explicitly Eq. (3.6) leads to 



hab,r = —2nXab + Ya\b + Yb\a 



(3.7) 



where Ya := habY''. Since expressions for the Ricci tensor components, Rij{;-y), of 7^^ 
are needed for some of the calculations we give these explicitly, in the present notation, 
in the appendix below. 

The Hamiltonian constraint (Eq. (2.9)) gives, using = 0, the following equation 
at dM: 

A -K. 



1,=,^^ (y)' 



(3.8) 



Since K > Q and we require O > for r < r_|_ this yields 



K 



rT 1'"='"+" ~"3 



(3.9) 



One is free to compute angular (but not radial) derivatives of such equations and 
deduce, for example, formulas such as 



I ^ J,a-U,... . 



(3.10) 



To get the next order Taylor coefficient of O, we compute the radial derivative of 
Eq. (2.9) and evaluate the result at dM to find 



n , 
— + 
n 



n 



n 
2' 



— —h°'^Xab 



lr=r-)- 



0. 



(3.11) 



The third radial derivative of Cl at dM can be computed by first taking the Laplacian 

of Eq. (2.9) and then reducing the resulting expression using the foregoing results. The 
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key formula resulting from this calculation can be written as 

1 



n 



r=r-|- 



1 



1 ~ 1 TT*"^ ^-J TT*"^ 

1 

+ 7; —{h'"'Xab),d 

2 n . ^^^^ 



(3.12) 



One needs the full, readily computed expression for the Laplacian of Q (given in the 
appendix) to evaluate the left hand side of Eq. (3.12) and, in particular, to extract the 
result for di^O, |^_^^ but, when only angular derivatives of this Laplacian are needed. 



one can use instead the limiting form 



(3.13) 



r=r+ 



Attempting to go beyond this level and compute d^VL 



r=r-|- 



however proves fruitless 



since, under further differentiation of the Hamiltonian constraint and evaluation at Scri, 
the coefficient of the fourth radial derivative of O at dM is found to vanish identically. 
On the other hand, given the results above for the first two radial derivatives of VL at 
the boundary, it is straightforward to evaluate the traceless part of the Hessian of Q, at 
dM and show that 



0, 



(3.14) 



r=r+ 



^0, 



r=r-(- 



and 



lr=r_|_ 



(3.15) 



K 



14 



r=r+ 



where, as usual, we write fl'^^ for Y^^^'^Q-im- Note that the quantity in square brackets 
in the final equation is just the traceless part of the second fundamental form Xai, 
induced on dM. 

By examining the angular components of Eq. (2.22) and comparing these with Eq. 
(3.15), we see that a necessary condition for regularity of the evolution equations at 
Scri will be the vanishing of the quantity 

li^a""^ := tt'""'' + ii^ih^'^h'"^ - ^ h°'^h'"^)Kd (3.16) 

on this boundary. This condition was derived by Andersson et al. in Ref. [12] and 
identified geometrically there as equivalent to the well-known requirement that the 
shear of X"*" should vanish. An examination of the remaining {ri) components of Eq. 
(2.22), together with a comparison of these to the first and second of Eqs. (3.14), 
shows that we shall need the additional regularity conditions tt*^^* = holding at Scri 
to avoid singularity there. We shall see momentarily however that these apparently 
new restrictions are in fact forced by satisfaction of the momentum constraints and 
thus do not represent additional limitations upon the free data. 

When the mean curvature is constant, as we have assumed, the momentum con- 
straint, Eq. (2.7), takes the form 



+ f L(7)7r'''^] - 2(a,f])7r*'-*^' = 0. (3.17) 



Evaluating this at 5M, and recalling that 0.^^ = leads immediately to the 

conclusion that 

TT*-^^* = (3.18) 

which, as noted above, are necessary for regularity at the boundary. Further important 

information results from computing the radial derivative of Eqs. (3.17) and evaluating 

the result at dM. The 'strong' form of the resulting equations (i.e., that attained before 
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the boundary conditions tt*^'^'"* = have been enforced) is given by 



tr ra 



2/^'"""n _tr rr , ^ab tr ab 



+ — 7T' 

n 



n 



n 



+ 2(^ + — A„5)7r*-« 
n n 



n n n J 



r=r-)- 



(3.19) 



and 



_tr CO -yc 
TT +n ( o-)\a AofoTT 



n . 



_2(!^ + ^^-^/,«^A,6)7r*-- 
n n 2 

+ 2(- n/i^'^A^b - -F^y'A„5 - — n,, + I^^Jtt*-^^" 

+ (- —n,r + h'"^Yd,r - —Y'^Y'^Xab 

n n 



4o. 



(3.20) 



r=r+ 



When the boundary conditions are imposed however, these yield simply 



I^trrry I A rAa6^tro6| | 

L ,r J \r=r+ I ^ J I 



r=r+ 



(3.21) 



and 



tr ca 



lr=r-(- 



(3.22) 



An application of L'Hospital's rule to the corresponding (apparently) singular terms 



n 



yields 



lim 

1 — >r-)- 



TT 



tr rr 



,7r 



tr rr 



n 



lini (-^) 

_ ^^A,,7r*'^-^ 



(3.23) 



K 
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I r=r_(_ 



and 



lim 

) — >r-)- 



TT 



tr rc 



^tr rc 

lim f — ^— ) 

\ „tr ab 



tr ca 



3 r .TT 



( ^2 )la 



(3.24) 



r=r+ 



Using the foregoing results for the derivatives of f2 we find, in the analogous way, reg- 
ular boundary expressions for the remaining (apparently) singular terms in Eq. (2.22): 



r=r+ 



hab 2 1 

{-^""lab-^ [^abXcdh'^'h''' + r?R''{^) - -i?(7) 

1 trlj trim -i i 



(3.25) 



6 (^^)2 
-^(^"^A.,),.]}|_^, 



2n 



r=r+ 



ad 



h 
n 



n\ab 



(3.26) 
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r=r+ 



and 



I — 

I r=r-)- 



^ n ^ 

^ n 

1 _ 1 1 ^trlj^trim yc 

- ^(/l"'A„6),.] +//-y^(/l^%d)[-/l^"/l^'Aab + /l"''Aab^^ (3.27) 

n 

„2 yehfa yfuea Y^Vf 

For the sake of generality, we have left Eqs. (3.25-3.27) in their 'strong' form - they 
could be slightly further simplified by imposing tt*^'"* = at dM. One could also 
substitute the explicit expressions for R^^ (7) given in the appendix, into the above, but 
since there are other Ricci tensor contributions to the full tt*"^*-^ evolution equations, it is 
more prudent to leave Eqs. (3.25-3.27) in their present form until the final evaluations 
are ready to be made. 

Equations (2.9) and (2.13) can be combined in an obvious way to yield the following 
equation for the conformal lapse function 

- ni^'m.ij) - Ni^'m.ij) - (3.28) 

+ + ^NQ -Jl^y^ lijlml 

= 0. 

Evaluating this at the boundary dM gives 



■ h°'^Xab 



Nr 1,=,/ {Y'N, - nN{'-^)} l^^^ (3.29) 
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where we have used Eq. (3.13) to reexpress the Laplacian of f2 at r = r+. Radially- 
differentiating Eq. (3.28), we can derive an expression for the Laplacian of N at 
dM, from which a formula for N^rr \r-r_^. ^® extracted. Using these results it is 
straightforward to evaluate the //^(V^V-'iV — |7*-'7"^"'VmVn-/V) contributions to the 
field equations at dM. One cannot compute higher radial derivatives of N in this 
way however since that would necessitate knowledge of the values of higher than third 
radial derivatives of Q at the boundary. In every case however, the Taylor expansion 
techniques allow the computation of precisely what is needed for evaluation of the 
evolution equations at Scri. We shall present these equations explicitly in the next 
section and show how they imply preservation of the regularity conditions tt^'^''* = 
and TT^"''' + n^ih^^W"^ - \h''^h^'^)\cd = at the conformal boundary. 

Though we have assumed constant mean curvature throughout, it would not be 
difficult to incorporate the contributions of variable mean curvature. For example, 
letting T = ^^7^ = , and allowing r to be variable, one gets for the momentum 

constraint 

QS/jTr^'"^^ - 2i}j7r''"'^ + \iij^''^diT = 0. (3.30) 
This leads to the modified regularity constraint 

- 07--9.r} 0. (3.31) 

Similar modifications are implied for tt^^'^ and for the radial derivatives of O and 

N at Scri. We do not, however, anticipate that such modifications would interfere 

significantly with the main conclusions derived herein. Instead, we believe that the 

CMC condition plays a rather inessential role in our analysis but we prefer to retain it 

because of the associated, partial decoupling of the constraints that it allows. 

We also note that, whereas our full constrained evolution proposal (as sketched in 

section II) entailed a constant scalar curvature normalization and harmonic coordinate 
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conditions for the metric 7jj, the corresponding (non-degenerate) eUiptic equations for 
the shift vector X^-^ (Eq. (2.15)) and the function F (Eq. (2.12)) played no role in the 
above analysis. Thus we do not believe that these particular choices are at all essential 
for the central conclusions derived herein. Rather we expect that many alternative 
formulations of the field equations could be adapted to putting the outer boundary at 
Scri and that our specific proposal is just one of many feasible possibilities for doing 
so. On the other hand, we also believe that many of our calculations have a rather 
'universal' character and will be applicable to a variety of alternative formulations which 
adopt different gauge conditions or normalizations from the ones we have chosen. 

We remark here that one of the main advantages of preserving, as we have, the strong 
forms of the relevant equations is that these lend themselves to further generalization 
(through modification of the gauge or normalization conditions or the introduction of 
material sources) by the straightforward computation of the additional terms necessi- 
tated by the desired modification. By contrast the information lost in passing to the 
weak forms of these equations would almost surely necessitate a rederivation of most 
of the relevant formulas from scratch. 

IV. Preservation of Regularity Conditions 

In this section, we fit together the various contributions to the tt*"^ evolution equa- 
tions derived above and show that they imply 

fl_.tr i^'i- f) 1 

^ S^""''"' + f^-y{h^'h'd - - h-'h'd)X,d) = (4.1) 

independently of any further restriction upon the geometrical data at Scri (i.e., without 

the need for any additional constraint on the boundary values of 

{hab, n, y*, N, Z", Xab, r or Rij{'y)}). We do this by first simply evaluating the 'strong' 

form of the relevant evolution equations and then noting, by inspection, that when the 
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regularity constraints, {tt^"^" = 0, fi^a""^ := tt^^""^ + fiy{h''''h^'^ - \h''^h^'^)\cd = 0}, are 
enforced these evolution equations reduce to the (trivial) forms given above. 

There is however, a subtlety in this seemingly straightforward procedure that we 
wish to address at the outset. This concerns the logical significance of the 'strong' 
forms of some of the relevant equations and the justification for replacing them with 
their corresponding 'weakened' forms in deriving the main results. To see the issue 
at hand, consider for a moment our derivation of the expressions (3.23 and 3.24) for 

tr vi 

the boundary values of the apparently singular terms using L'Hospital's rule. We 
chose, for simplicity, to present those results in their 'weak' forms by dropping all those 
contributions from equations (3.19 and 3.20) that vanish at Scri by virtue of the reg- 
ularity conditions even though the retention of such terms would still have yielded a 
completely regular result. The strong forms of such evaluations, though perfectly regu- 
lar at Scri, entail a certain logical contradiction - one assumes the regularity conditions 
hold in order to extract finite limits from otherwise singular expressions but then, after 
applying L'Hospital's rule, leaves in the tt*^*"* and a"-^ contributions as though the reg- 
ularity constraints were actually being relaxed. For other purely regular terms in the 
evolution equations, no such contradictory assumptions are involved in their evaluation 
and so their strong forms seem less problematic. 

Of course, one is always free to drop all weakly vanishing terms and restore logical 
consistency so no harm is done in retaining such terms in early stages of the calculation. 
We suspect however that in some eventual, deeper mathematical study of these issues, 
it may be important to know the actual structure of such weakly vanishing 'forcing' 
terms and so we have uniformly retained them in the derivations reported below. 

The reader may well wonder, however, whether we are not then obligated to prove 

a strong version of conservation of the regularity conditions by establishing a form of 

hyperbolicity of the boundary evolution equations or (essentially equivalently) deriving 
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suitable 'energy' estimates to show that these equations have only the trivial solution 
for vanishing initial data. Here, however, the illogic in the derivation of such strong 
equations comes to the foreground. If one imagines that the regularity conditions 
are not necessarily enforced at Scri (as would be implicit in the use of the strong 
form of the evolution equations) then one has no logical right to make contradictory 
assumptions for the 'evaluation' (using L'Hospital's rule) of genuinely singular limits 
such as lim and lim On the other hand, the explicit contributions of these 

seemingly problematic expressions to the corresponding evolution equations, namely 
the terms 

take the forms of 'frictional forcing' terms that (having good sign for evolution towards 
the future) seek to drive the quantities tt*"^^* and a"''' towards vanishing values with 
effective frictional coefficients that diverge like as Scri is approached. Intuitively 
speaking then, the regularity variables {tt^'^'^*, a""^} feel precisely zero force when they lie 
on the 'constraint submanifold' defined by the regularity conditions Tr**^*^* = 0, cr'*^ = 
but would be subject to an infinitely strong 'frictional' restoring force if they tried to 
'float off' this manifold. Finding a precise mathematical characterization of this state 
of affairs might prove to be a subtle task, but we are encouraged that the indicators 
that we have identified all point towards stable preservation of the regularity conditions 
at Scri. 

In the previous section, we showed how to evaluate the most problematic contribu- 
tions to the evolution equations at Scri. The remaining terms (i.e., those involving the 
Lie derivative of tt*"^ , non-singular algebraic expressions and terms involving the Ricci 

tensor of ^ij) are straightforward evaluations that require little comment. We recall 

: d 



however, from section II, that the angular components of the shift vector, X^t^j, have 
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the following form at Scri 



I = {—N-f'^m i + Z'^]\ (4.2) 

lr=r-|- ^ fC lr=r-|- ^ ' 



or equivalently, using foregoing results on the behavior of fl at dM, 



I =1— y^ + Z^jl . (4.3) 



Accordingly, Lie derivatives with respect to the 2-dimensional vector fields 

(2)^ := (2)^ := za_|_ and (2)y := are expressed as C,2,x,jCi2,z and 

>C(2);P respectively. 

Assembling the various elements of the (strong forms of the) tt*"^^* evolution equa- 
tions, one thence obtains 



1^^*'''} \r=r+^ {-X^y-^- - 2Xiy-'' + (XV^^^),a (4.4) 

2N 

[(tt j (n + YcY ) + 2Yan tt + habT^ tt J 



"I Fa "I '^afeTT +( 1 Y Y Aab + nil AabjTr 

n ' n n n n 

" 23 



2/V 

r^tr rr _tr r(i/„2 , "5^ xac\ , _tr rr _tr da-iA 
[TT TT + y j + TT TT la 

+ TT TT ra + /labTT TT J 

+ ^ [_ 2(^ + ^ - Ih'^^Xa.y^"^ 

n n ' ' 

^ n n 

~ , .TT*''''' + ii^(h'"^h'"' - lh'"^h''^)\ab. 
+ [ -z ]|c 

Proceeding in the same way, one finds for the evolution equation of the densitized shear, 
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|{(7r*-^ +/.,(/^-/^^^ - h'^')Xa,)} (4.6) 



A r-2N 



^ac \\c 1 / h^'^Xi 



^,iV(F^/^'^/ + y//^-)[(^)'^-±(:-^)J 



n ' 2 n 
/x^ 2n ^ jiy fly I** 

2iV 

+ 27V[^7r*-/(/.'^''Aa6) - \x'fhab^''"'' + //^(2A^«A{ - ^A'^^/i«^Aa6 - \h''^ X""' Kh)] 

+ -^(7r^-^ + /.^(/.^«/.^''- J/.«^/.'^'')A„5) 
n ^ 2 ' 
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For the first two equations (4.4) and (4.5) it is clear by inspection that the right 
hand sides vanish 'weakly' (i.e., vanish when the regularity constraints tt*^^* = and 
fx^a''^ := [7r^'''^+ii^{h'"'h'"^-^h''^h'"^)Xcd] = are enforced at dM). To verify that the 
right hand side of the third equation (4.6) also vanishes weakly, requires a little work. 
For this purpose it is important to note that since the 3-dimensional trace, 7ij7r*^*-', of 
^trij vanishes identically, one has 

7,^.7r*-^- = habTr"""' + (n^ + y,y'=)7r*'-^- + 2y«7r^'^^« (4.7) 
= 

and thus finds that the 2-dimensional trace, /labTi"*"^"^) of n^^"'^ vanishes weakly at dM. 
It is also useful to recall that an arbitrary 2-dimensional, traceless, symmetric tensor 
s°'^ satisfies the identity 

s'^'^s^'^Kd = lh''\h,ehdfs'"'s^f) (4.8) 

when s"'' = s''" and habs''^ = 0. 

Thus equations (4.4 - 4.6) reduce weakly to equations (4.1) and so imply the preser- 
vation of the regularity conditions throughout the evolution at Scri. It is important to 
note, especially if alternative gauge conditions, normalizations or formulations of the 
field equations are under consideration, that our results involve no implicit restriction 
upon the boundary values of the geometric data {hah, F", A^, Z", Aab, T or Rij{j)}. 
Our main tools have been simply the enforcement of the Hamiltonian and momentum 
constraints together with the straightforward application of Taylor expansions and 
L'Hospital's rule for the evaluation of apparently singular terms at dM. Although we 
have also imposed the CMC slicing condition (and peripherally, the spatial harmonic 
gauge condition), we do not believe that these were at all essential for our principal 

conclusions and that one could relax them without harmful effect. 
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V. Alternative Evaluation of Singular Terms 



As is well-known [14, 15] the electric components of the Weyl tensor, for a solution 
of the vacuum field equations, can be expressed in terms of the physical Cauchy data 
{9ij ) ^ij ) ^ 

Eij = Rij{g) - K,mK^^ + KijiTgK (5.1) 

where K"^- = g^^Kij and tigK — g^^K^i. Looking at this equation though, one may 
well wonder why, since the Weyl tensor is conformally invariant, the above formula for 
Eij is not. The answer in part is that the above expression does not correspond to the 
fundamental formula for Eij but rather to a representation of it that has been trans- 
formed, through the application of the ADM field equations, to a form in which time 
derivatives have been eliminated in favor of spatial ones. Thus the above expression 
for Eij has, in effect, inherited the failure of the ADM equations to be conformally 
invariant. But this same lack of conformal invariance in the ADM equations is pre- 
cisely the feature which led, upon conformal rescaling, to the appearance of the singular 
terms that we have been concerned with. Thus it should perhaps not be surprising to 
find that there is a close relationship between those singular terms and the conformal 
transformation properties of the above expression for E^j . 
To see this explicitly, let us define 

■■= f^gEr^ = ^^99''Eml (5.2) 

and reexpress the physical variables in terms of the conformal ones introduced in section 

II. Even though the trace of £^ vanishes by virtue of the Hamiltonian constraint, it is 

convenient to write the resulting formula in terms of the explicitly trace-free quantity 

~ \^rn^k sincc this facilitates comparison with our earlier derivations. Using the 

well-known conformal transformation properties of the Ricci tensor, and introducing 
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the gravitational momentum variables tt^'^*-' in favor of Kjj, one easily arrives at: 
_ 2 KN _ ~ V^W^Q _ 1 y^V^. 

3 Q n z n > 

= 2iV^^(i?^-'(7)-l7^-'i?(7)) 

The left hand side of this equation consists of precisely the apparently singular terms 
in the tt*^*-' evolution equation (2.22) whereas the right hand side consists of purely 
regular terms provided that is regular at T+ . 

The reader may wonder however, whether we have hidden some singular behavior in 
the notation by choosing the mixed, densitized form S^^ to represent the electric compo- 
nents of the Weyl tensor. That we have not done so however, may be seen by writing 
down the evolution and constraint equations satisfied by £^ and the corresponding 
magnetic components B^, expressible in terms of ADM variables as 

BJ, := ^,B^ = ^igg^'Bmi (5.4) 
= ^ e'^^{Vi{g)Kim-Vi{g)Kim) 

or, upon assuming that tigK — constant as above and reexpressing K^j in terms of the 
gravitational momentum, as 

Brii^ — 9mn9is{'Vl{g)7T''''^). (5.5) 

These Maxwell- like evolution and constraint equations, derivable directly from the 

defining formulas for Sj^ and Bj^ through an application of the ADM equations, can 
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(5.3) 



be expressed as: 

d_ 
di 



S/ = (C^S)/ (5.6) 



'^rsi ■'■^ ,m I '^j 
5 N-fjk p / trfcm 

+ 2 ^^-^ 

N'Jkm trfcnc m 



1^1 

1 N 



I ^ " ' ^, ^tr kl c m 

2 ^ 



2 //^ 



V,(7)^r^' = - Smlr^'^'^ B/, (5.8) 



^M^J. = + ^mlr7r"''' S/, (5.9) 

where X := X'^-£^ and signifies Lie differentiation witfi respect to X. 

Tliougli we liave written tlie above in terms of tlie conformal metric variables, we 

could just as well have replaced them by the physical metric components since it is 
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easy to see by inspection that all of Eqns. (5.6-5.9) are conformally invariant. For this 
reason, therefore, the possibility that some regular null hypersurface in the unphysical 
geometry may in fact play the role of future null infinity in the physical geometry is 
completely invisible to the quantities {£^,B^), which evolve without any interaction 
with the conformal factor 0,. Thus, unless some singular behavior is put 'by hand' 
into the initial conditions for and Bj^ at X+, there is no reason to anticipate 
that these quantities will blow up along null infinity, at least for sufficiently short 
time intervals during the evolution. This plausibility argument that singular boundary 
behavior for {S^,Bj^) should not be put in 'by hand' can, however, be strengthened to 
a mathematical proof (due originally to Penrose [2]) that sufficiently smooth conformal 
compactifications must necessarily have vanishing unphysical Weyl tensors at Scri. As 
we shall see below, this result in turn implies that ^Sj^ and ^BJ^ must both vanish at 
Scri. 

Thus the issue of whether ^Sj^ should be assumed to vanish at Scri hinges ulti- 
mately upon the 'reasonableness' of Penrose's smoothness hypotheses (which require 
C^-differentiability of the conformal metric out to the boundary). But this question 
opens the Pandora's box of concerns as to whether one should allow so-called 'poly- 
homogeneous' (or polylogarithmic) boundary behavior as a 'natural' alternative to 
Penrose's smoothness requirements. The issue of whether polyhomogeneous Scris are 
actually needed for sufficient physical generality is still an open question upon which 
we currently do not wish to express a definite opinion. For simplicity therefore, we 
have adopted the traditional Penrose viewpoint that null infinity should be sufficiently 
smooth that the corresponding argument for the vanishing of the unphysical Weyl ten- 
sor at Scri can be applied. As we shall discuss more fully below however, most of our 
calculations are actually compatible with the possibility that Scri could be polyhomo- 
geneous instead of smooth in the Penrose sense. 
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Assuming therefore that the quantities fiSj^ vanish at T+ we obtain from Eq. (5.3) 
an independent expression for the apparently singular terms in the evolution equations 
for TT*^ '-^ in terms of regular quantities. How do these compare with our previous calcu- 
lations? It is straightforward to verify that, for the ri components of these expressions, 
the two formulas agree weakly (i.e., upon satisfaction of the regularity conditions) at 
the conformal boundary but not strongly. A direct comparison of the expressions for 
the angular components however, is not possible since our previous calculation, through 
its application of L'Hospital's rule, required the boundary values of ir^^,^^ whereas the 
Weyl formula given above involves no derivatives of tt*"^"^ at X"*". L'Hospital's rule was 
of course also used for the evaluation of the ri components but there we could use 
the momentum constraints to compute the corresponding radial derivatives, tt*^/"* at 

. We did not need boundary expressions for the quantities Tr*^^,"^ in the calculation 
to show that the vanishing of the shear of Scri is preserved since the contributions 
of these radial derivatives (including the Lie derivative of tt*^^"'' and the time deriva- 
tive of Xab) actually conspired to cancel one another. But this earlier calculation only 
succeeded to give an expression for the time derivative of the shear at Scri and did 
not yield a formula for either TT^y'' or {^^{h'^'h'"^ - \h''V)\cd),t separately at the 
conformal boundary. To obtain the latter we require some new input. To derive an 
expression for tt^^^ at Scri, and thereby to complete the derivation of a formula for 
TT*^^"'' at this boundary, we equate the two expressions for the angular components of 
the singular terms at Scri (given by Eq. (3.27) and by the angular components of Eq. 
(5.3)), assuming that they must agree, at least weakly, for any regular solution of the 
field equations. The resulting formula is: 
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n 



Trtr ef tr ef 



(5.10) 



(Aa6-i/ia6(/i=''Acd)) 



n 16 



lx+ 



^0 



where the (weak) vanishing of tt*"^^' \ and hat'^^^"'^ 1 7-+ has been imposed to simphfy 



the expression. 

An alternative derivation of the above formula can be developed by considering the 
conformal transformation properties of the ADM expression (5.5) for the magnetic 
components of the Weyl tensor in the analogous way. Reexpressing this formula in 
terms of conformal variables, one arrives at 



(5.11) 



,7r 



tr rs 



= e'^'l,r V,(7)(^) 

The antisymmetric projection, Sksiil^'^^Bj^), of the above expression vanishes by 
virtue of its equivalence to the momentum constraint so one needs only to consider 
the symmetric projection V*' defined by 



2V' 



(5.12) 



TT 



tr ks 



e'^'lik V,-(7)(-— ) 



+ e^^'^7.fe V,(7)(^). 
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Evaluating the right hand side of this formula and appealing to the regularity con- 
straints TT*"^" = and [Tr*-^"^ + - ^h'"'' h"'^) Xcd] = as well as to Eq. (5.10) 
and the momentum constraints for the computations of 7r|^"^ \x+, one can show by a 
straightforward calculation that 

=0 (5.13) 

as would be expected from the assumption that fiB;™' vanishes at Scri. In other words, 
the vanishing of at Scri necessitates the enforcement of Eq. (5.10) at this bound- 

ary. 

Finally though, we can remove the guesswork in the aforementioned plausibility 
argument by appealing to Penrose's well-known result that states that the unphysical 
Weyl tensor, W"p^g{^'^^^), should vanish at Scri [2,3]. Translated into our notation this 
result gives directly that ft £^ and Q, BJ:^ vanish there. To see this, recall that by 
conformal invariance, 

and thus that 

W^^,s{^^^g) = n,,((^)^) (5.15) 

Letting n^-^^ be the (physical) timelike unit normal field to the chosen slicing (so that 
^^^g^^n'^'ri'^ = —1) we get 

W^p^si^^^ 9)71(^71' (5.16) 
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where n^^ = ^ yields the corresponding unit normal field defined relative to the con- 
formal metric (so that '■^^7^t/n'^ n'^ = —1). 

The electric components, Eij, of the physical Weyl tensor are given (in coordinates 
for which the chosen slices coincide with = t = constant hypersurfaces) by 



Eij = W.^jsi^^^g) n' 



(5.17) 



and thus vanish at Scri by the aforementioned argument. Since, however. 



mj 



(5.18) 



= 7 



mj 



it thus follows that flSji^ vanishes at Scri. By taking appropriate duals of the above 
expressions it follows in the same way that Bj^ also vanishes there. Either one of 
these results, as we have seen, is sufficient to imply Eq. (5.10) for the radial derivative 
of TT^'^"'' at the outer boundary. 

Though somewhat peripheral to the above discussion, we conclude this section with 
some remarks on the evolution of geometric data along the conformal boundary and on 
the choice of a "conformal gauge" at Scri. Recalling that the shift vector components 
X' satisfy (c.f., Eqs. (2.18), (3.9), (7.4)) 



one sees that 



I r=r^ 



N 



n 



r=r+ 



lr=r+ ] n 



(X^ yc ^ j^c^ 



r=r+ 
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(5.19) 



r=r-(- 



(5.20) 



r=r+ 



Using this, it follows from a direct evaluation of the angular component of ^ij^t at 
Scri that the metric hah induced on 2-dimensional slices of Scri satisfies the evolution 
equation 



b,t 



tr cd 



(5.21) 



2/V 



+ Z'^hab,c + Z'^a^cb + Z'^b^ac} 

which clearly is weakly equivalent to 



r=r-(- 



hab,. 



r=r+ 



AT" "p 

= {y M^ + 3(/.^^Ae/)) 



(5.22) 



r=r+ 



One could thus exploit the freedom to choose Dirichlet data for the function F (which 
heretofore has remained unconstrained) to arrange that (r+3iV h'^^ Ag/) |^_^^ vanishes 
on each t = constant slice of Scri. This choice would reduce the above evolution equation 
to the essentially trivial form 



ab,t 



= {-^(2)^ hab} 



r=r-|- 



(5.23) 



r=r+ 



or, if one also exploits the freedom to set 
form 



= at Scri, to the manifestly trivial 



r—r^ 



h 



ab,t 



(5.24) 



r=r-(- 



when Z*^ 



= (r + 3iV h'^fXef) 



r=r+ 



^0. 



r=rj^ 
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r=r_|_ 



A 'covariant interpretation' of this (potential) choice for (F + 3N h'^^ X^f] 
be uncovered by applying the (conformal) wave operator 0(4)^ to the conformal factor 
fl and evaluating the result at Scri. A straightforward calculation yields the result 



1+ 



A r , 1 r 



(5.25) 



r=r-|- 



Thus the special choice of Dirichlet data for F discussed above corresponds to imposing 
the 'conformal gauge condition' that 0(4)^ ^ vanish at Scri. By appealing to the 
conformal transformation properties of the 4-dimensional Einstein tensor, one can show 
that this conformal gauge choice further implies that 



^0 



(5.26) 



i.e., that the full 4-dimensional conformal Hessian of Q vanishes at Scri [16]. Since this 
result simplifies numerous formulas evaluated at Scri, it is often imposed to streamline 
various calculations at the conformal boundary. The above discussion shows that such 
a choice of conformal gauge is fully compatible with our formalism and determines a 
correspondingly 'natural' choice for the Dirichlet data for F while leaving unconstrained 
the boundary values of N and Z'^. 

Assuming Penrose's result that the quantities and should vanish at Scri 

we can use Eq. (5.10) to eliminate the radial derivative, tt^^^'^'^ from the equation 
of motion for tt*^^-^ \^_^_. In view of Eqs. (4.6) and (5.21) however, it clearly suffices to 



present instead the (slightly simpler) equation of motion for Xab A straightforward 
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calculation making use of Eqs. (2.21), (3.2) and (3.7) yields 

{2n^} |_^^ (5.27) 

fij 6 
fij fly 

+ [ZYa — n TT H n hac^ \\b 

-{ — /^ac/^6d7r^'^^' + — [l^an^r*'^'-^ 

- 2NXab,r + 2N{C^^)y\ab) + 2n{C(^)z^ab) 

2N 2N 

+ r«/i6d7r*'^^'' + n/iad7r*'^'^')] 

- [^'(^)|b]|a- [n'(^)|a]|b 
+ i/.«5(l^Te-r,)} 1^^^^. 

For completeness we also list the equations of motion for the metric functions n and 

yc. 

n at 3 



r? 2 

+ 2X'=^ - 2r^x': + -(nxn ^ , 



^ = 2^^— nV*'^'"^ + —n\''^^ (5.29) 

+ x"y;^ - x^'^^r'^ + x;„(n2/t"^ - r"r=) 
+ (x'"y'^ + x^),^. 

Notice though that even when these are evaluated at T+ they involve radial derivatives 

of the unknowns as well as those of X* at Scri. The same is true of Eq. (5.27). Thus, 
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unlike Eqs. (4.1) or (5.21) for example, Eqs. (5.27) - (5.29), when evaluated at Scri, are 
not fully intrinsic to this boundary. This would of course be all the more true if Pen- 
rose's result is not presupposed and the radial derivatives, n^r'^-^ \J_^., can therefore not 
be eliminated as discussed above. On the other hand, if Penrose's smoothness assump- 
tions are made, then one can return to the momentum constraints (3.17) and compute 
an additional radial derivative thereof in order to derive expressions for tt*^^*^ |^_,_ and 
TT^rr"' lx+" Thesc may be useful in extending the Taylor expansions for the quantities 
^trrt about T+ but one should keep in mind that the additional smoothness require- 
ments needed to justify this calculation may be incompatible with the generality desired 
for the solutions under study. We shall return briefly to this issue in the concluding 
section. 



VI. Concluding Remarks 



In the section above we mentioned the (still controversial) issue of deciding whether 
Penrose's smoothness hypotheses are overly restrictive from the stand-point of physical 
generality. Fortunately however most of our central results do not require stringent dif- 
ferentiability assumptions and are in fact compatible with so-called 'polyhoniogcneous' 
(or 'polylogarithmic') behavior at Scri. To see this we appeal to the conclusions derived 
in Refs. [12,13] wherein one finds that polylogarithmic behavior, if present, sets in at 
the level of the fourth radial derivative of fl and the second radial derivative of tt*^*-' 
at . But our Taylor expansion calculations only determined the radial derivatives 
of fl up to the third order and those of tt*^''* up to the first. Only if Penrose-level 
regularity was additionally assumed could one go further and compute the quantities 
{tt*^^'', TT*^^", TT*^'*''} 1 for example, but none of these were needed in the demonstra- 
tion of regularity at the boundary or its preservation under time evolution. 

If polylogarithmic Scris are allowed, then one must forego the argument that the 
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unphysical Weyl tensor vanishes at null infinity and therefore also the corresponding 
evaluation of the aforementioned, higher radial derivatives at since the latter, in 
general, may not have well-defined limits at this boundary. But since none of these 
were essential to our central results, we believe the latter are fully compatible with 
poly logarithmic, as opposed to smooth, behavior at Scri. 

Some numerical relativists may wish to put their outer boundaries at Scri but prefer 
not to adopt constrained evolution or prefer to avoid the use of our CMCSH gauge 
conditions and wonder, accordingly, whether our results have any relevance for them. 
While we cannot draw definitive conclusions that are sure to apply to an arbitrary 
numerical setup we can nevertheless make several general points in this regard. First 
of all, the main tool needed for verifying finiteness and then actually evaluating the 
apparently singular terms at Scri has been the strict enforcement, at least to the requi- 
site order in a Taylor expansion, of the Einstein constraint equations and the vanishing 
shear condition at the outer boundary. Thus we would anticipate that any strategy 
for say free (as opposed to constrained) evolution would have to be coupled with strict 
enforcement, to the needed order, of these constraints at the boundary. Otherwise, 
there seems to be no hope for showing that the evolution equations are actually regular 
at Scri. On the other hand, the 'universal' character of our calculations (i.e., the fact 
that they hinged purely on imposition of these constraints) would seem to show that 
they apply equally well to a variety of formulations of the field equations and not just 
our own. 

Furthermore, gauge conditions played a rather peripheral role in our analysis and, 

for that reason, we do not think that our particular choice was at all crucial for the 

main results. We have already sketched how the CMC condition could be relaxed and 

our calculations modified accordingly. If the mean curvature, tr^X = gijK'^^ , is allowed 
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to evolve, its equation of motion, in our notation, can be written 



^(tr,K) = -^{Y^N,ij) - N{^^m.,,^) - ^QN (6.1) 



3 -77^""^^^'^^ 

-Nn — - — 

+ X^(tr,K),, 



and has no irregular behavior at Scri. Furthermore, our preservation of regularity 
calculations did not make use of our elliptic equation for the shift vector X^-£^ (imposed 
to preserve spatial harmonic coordinate conditions for jij) or of our normalization 
condition R{^) — constant. In addition, they did not lead to any restriction upon the 
asymptotic gauge data {N, Z"'-^} \^_^_. They did make use of our elliptic equation for 
the conformal lapse function N but this was forced in large measure by our imposition 
of CMC slicing conditions and use of a simplified form of the momentum constraint 
which only holds in CMC gauge. We have little doubt that the preservation of regularity 
calculations would work equally well in a variety of other frameworks. 

Concerning the issue of a stable, numerical implementation of the regularity condi- 
tions at Scri, this has not yet been carried out. A key point though is that, in any such 
implementation, the computer would not be expected to re-derive (by delicate limit- 
ing procedures) regular expressions for the apparently singular terms in the evolution 
equations at the boundary. These would be explicitly provided as boundary conditions 
through formulas such as (3.25)-(3.27). In view of the comments in section IV regard- 
ing the 'frictional restoring forces' in effect near Scri, it is very plausible that a stable 
implementation of these boundary conditions is numerically achievable. 

In summary, we see little reason why a wide variety of numerical setups could not 

consistently place their outer boundaries at Scri. As a prominent numerical relativist 

recommended to one of us some years ago "If you want to go to Scri - just do it". 
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VII. Appendix 



In the notation we have used, the Einstein tensor, 



G^^(7):=^^^(7)-^7^^'%), (7.1) 



of the metric 7^ has components 



+ ^^^R{h)l (7.2) 
+ ^ [XahXcdh'^'h}"' - {h'^'^Kd? + ^^^R{h)\ , 

^ 17 j = Aafe,r - i-'(2)Y^ah 

\ n n J 

- n\ab + n ^^^Rab{h) + 2nh'"^XadXbc - nXabih^"^ Xcd)] 

+ \{h''^ - ^^) [XahXcdh'^'^h'"' - {h'^'^X^af 
+ (2)i?(/i)] 



where ^'^^Rab{h) and ^'^^R{h) are the Ricci tensor and scalar curvature of the 2-metric 

habi I a signifies covariant differentiation with respect to this metric and >C(2)y signifies 

Lie differentiation with respect to the 2-dimensional 'shift' field ^"gfs: = ^"^^Y . The 
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Laplacian of a function F is given by 



7^^F;,,- = {h-^ + ^) - ^{h-'^Xab) (7.3) 



2 

- -Y'^F,,h''Xab + -^F,Y]^ 
n ' I" 

TP VAO 21^'* 

IF F 

~l~ oF,rr ~l~ '^.r o"^ ^^,r 



The contravariant components, 7*-' , of this metric are given by 



1 y 

= ^, 7'"" = - — , (7.4) 

^ab ^ ^ab ^ 



(where h is the contravariant form of hab) and the volume element is expressible 
as 

jj.^ = nfih (7.5) 

where /ih is the area element of hab- 
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